The Incompressible Navier-Stokes Equations From Black Hole Membrane Dynamics 
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We consider the dynamics of a d + 1 space-time dimensional membrane defined by the event 
horizon of a black brane in (d + 2) -dimensional asymptotically Anti-de-Sitter space-time and show 
that it is described by the d— dimensional incompressible Navier-Stokes equations of non-relativistic 
fluids. The fluid velocity corresponds to the normal to the horizon while the rate of change in the 
fluid energy is equal to minus the rate of change in the horizon cross-sectional area. The analysis is 
performed in the Membrane Paradigm approach to black holes and it holds for a general non-singular 
null hypersurface, provided a large scale hydrodynamic limit exists. Thus we find, for instance, that 
the dynamics of the Rindler acceleration horizon is also described by the incompressible Navier- 
Stokes equations. The result resembles the relation between the Burgers and KPZ equations and 
we discuss its implications. 

PACS numbers: 11.25.Hf, 47.10.ad, 11.25.Tq 
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The dynamics of fluids has remained an unsolved prob- 
lem for centuries. Unlocking its main features, chaos and 
turbulence, is likely to provide an understanding of the 
principles and non-linear dynamics of a large class of sys- 
tems far from equilibrium. The fundamental formulation 
of the problem of the nonlinear dynamics of fluids is given 
by the incompressible Navier-Stokes (NS) equations 

dtVi + v~djVi = —diP + vdjjVi + fa , 



given 
1,1 

(1) 



where Vi(x,t),i = l,...,d, (d> 2) obeying diVi = is the 
velocity vector field, P(x, t) is the fluid pressure divided 
by the density, v is the (kinematic) viscosity and fi{x,i) 
are the components of an externally applied force. The 
equations can be studied mathematically in any space 
dimensionality d, with two and three space dimensions 
having an experimental realization. Among the central 
questions posed by the NS equations is the existence of 
singularities in the solutions and the statistics of the so- 
lutions in the limit of small v, both with and without 
forcing. 

The scarce progress achieved in understanding the NS 
equations prompts one to look for other frameworks of 
viewing the dynamics that might offer new insights. In 
this letter we explore a new viewpoint that can lead to 
a development of geometrical methods to study the NS 
equations. We represent the whole spatio-temporal pic- 
ture of the velocity field in terms of a structure of a null 
hypersurface (whose normal vector is also a tangent vec- 
tor), embedded in a bulk space-time and evolving accord- 
ing to the General Relativity field equations. Thus, the 
dynamics of the NS equations is related to the dynamics 
of the geometry, as described by the relativistic Einstein 
equations. 

A particular case of a geometric representation of the 
NS equations has been emerging recently based on the 
AdS/CFT correspondence relating conformal field the- 
ories (CFTs) in a flat (d + l)-dimensional space-time 
to gravity (string) theory on asymptotically (d + 2)- 
dimensional Anti de Sitter (AdS) space-time [3| (for a 



review see Q). Here the spatio-temporal pattern of fluid 
motion is mapped onto a black brane solution (a black 
hole with planar topology) in an asymptotically AdS 
space-time. The map is based on the realization that 
at large spatial and temporal scales the CFT dynamics 
reduces to hydrodynamics, and the dual gravity descrip- 
tion consist of long wavelength, long time perturbations 
of the black brane solution [f|. 

This correspondence has become relevant to the study 
of ordinary fluids since it has been shown that hydro- 
dynamics of any CFT, although relativistic intrinsically, 
does admit a non-relativistic (slow motion) limit given by 
the incompressible Navier-Stokes (NS) equations describ- 
ing flows of non-relativistic fluids [s |7|, l8[ . Thus, applying 
the correspondence between the black brane solutions of 
gravity and the flows of the CFT in the non-relativistic 
sector of the latter, allows one to construct a geometrical 
description of the NS equations. 

The way the black brane horizon geometry encodes 
the boundary fluid dynamics is reminiscent of the Mem- 
brane Paradigm in classical general relativity, according 
to which any black hole has a fictitious fluid living on its 
horizon [9[. In the AdS/CFT correspondence, the real 
fluid whose dynamics we wish to study is at the bound- 
ary of the space-time. It is natural to ask to what extent 
it can be identified under the duality map with the mem- 
brane paradigm fluid (see e.g. fiol. Illl|). 

We will consider this question in a rather general setup. 
Us ing, the Membrane Paradigm approach as developed 
in [l2j we will first analyze the dynamics of a membrane 
defined by the event horizon of a black brane in {d + 
2)-dimensional asymptotically AdS space-time, showing 
that it is described by the incompressible NS equations 
of non-relativistic fluids. The fluid velocity corresponds 
to the normal to the deformed event horizon, while the 
rate of fluid energy change is minus the rate of change in 
the horizon cross-sectional area. 

The analysis that we will perform also holds for any 
non-singular null hypersurface when a large scale hydro- 
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dynamic limit exists. Thus we will find, for instance, that 
the dynamics of the Rindler acceleration horizon is also 
described by the incompressible Navier-Stokes equations. 

The connection between the horizon hypersurface dy- 
namics and the NS equation is analogous to the con- 
nection between the Burgers [l3| and the Kardar-Parisi- 
Zhang (KPZ) equations [14j . Our result shows that real 
turbulence may also be seen as resulting from a physically 
natural surface dynamics. 

In the following we will use the convention 8irG — c = 
% = fee = 1- We will consider a (d+ 2)-dimensional bulk 
space-time M with coordinates X A , A = 0, d+ 1 with 
a Lorentzian metric gAB- Let H be a (d+ l)-dimensional 
null hypersurface (notion akin to horizon) characterized 
by the null normal vector n which components n A obey 

n ■ n = g A Bn A n B = . (2) 

Note that this condition implies that for a null hy- 
persurface the normal vector is also a tangent vector. 
We define the hypersurface in the bulk space-time by 
x d+1 = r — const, and denote the special coordinate sys- 
tem of [IH as x^ = (i,a; l ),z = 1,...,<± The coordinate 
t parameterizes a slicing of space-time by spatial hyper- 
surfaces and x 1 are coordinates on sections of the horizon 
with constant t. In this coordinate system n r = and one 
may choose the normalization n* = 1. We identify the 
remaining components with the fluid velocity in Eq. |T]) 

n l = v l . (3) 

In the case of black branes in AdS, the event horizon 
of a static, unperturbed, equilibrium solution is located 
in the bulk space-time at r — ttTq, where To corresponds 
to the Hawking temperature. The horizon coordinates 
(t, a; 4 ) can be identified with time and space Eddington- 
Finkelstein (EF) coordinates in the AdS boundary. We 
consider slowly varying (long wavelength, long time) per- 
turbations of the static black brane solution, so we work 
to leading order in an expansion in horizon derivatives. 
We also consider the slow motion limit where v l is a small 
perturbation. Temporarily restoring c, the horizon coor- 
dinates become (ci, x l ) and the non-relativistic slow mo- 
tion limit corresponds to v 1 /c -C 1. In what follows, in 
order to keep track of the different terms we impose the 
scaling d t ~ e 2 , v 1 ~ di ~ e, where e is a small parame- 
ter corresponding to c" 1 [22|. To 0(e) the (d + 2) black 
brane metric in EF coordinates is 

d 

ds 2 = -r 2 fdt 2 + 2dtdr + r 2 dx l dx % 

t=i 

2tt 4 T 4 

^-v t dx l dt - 2 Vi dx l dr, (4) 

ir 4 T 4 

where / = 1 pV 2 -. 

The first fundamental form of the horizon is the space- 
time metric restricted to it. In the horizon coordinate 
system it takes the general form 

ds 2 H = hi^dx 1 - v l dt){dx j - v j dt) , (5) 



where h%j is the metric on sections S of the horizon H 
at constant t. In our case, we work to leading order in e, 
considering only the terms zeroth order and linear in v l . 
For the black brane, at zeroth order 

h[f = (7tT )% . (6) 

The details of the subleading term, which is of order e 2 , 
will not be needed for our analysis. 

The second fundamental form of the horizon hyper- 
surface is the extrinsic curvature. In order to construct 
it consider the space-time covariant derivative Va pro- 
jected into the horizon surface and acting on the normal 
vector n A . Since n ■ n = ^ we have in our coordinates 

n A y„n A = . (7) 

This implies that V M n^ is tangent to the horizon and 
can be expanded in the horizon basis e A 

V,n A = K^e A , (8) 

where K" acts as the extrinsic curvature of the horizon. 
Together, the first and second fundamental forms pro- 
vide a complete description the embedding of the null 
hypersurface in the bulk space-time. 

Consider Lie transport of along the null normal 
vector n, which is given by the Lie derivative C n 

Cnhij = idthij + ±(^T ) 2 (A^ + DjVi) , (9) 

where Di is the covariant derivative with respect to the 
metric Note that the i,j indices are lowered and 

raised by the horizon metric hij in The above ex- 
pression can be split into its trace part (the expansion 0) 
and trace free part (the shear o~ij) 

= ^h'uhh,, ■ D.r' , (10) 
= C n hij - 9hij / d . (11) 

For the black brane we get to leading order the 0(e 2 ) 
expressions 

= dm , (12) 

°w = k^To) 2 (diVj + djVi - 2d k v k 5 i] ld) . (13) 

Using n and e A as a tangent basis, the components of 
the horizon extrinsic curvature are generally 

K = k(x) , (14) 
K? = *h , (15) 
K) = o) + 05)/d. (16) 

k(x) is the surface gravity defined by n B V B nA 
■ ." > We parameterize it as 

k(x) = 2ttTo(1 + P(x) ~ v 2 /2) , (17) 

where P(x) and v 2 scale as e 2 . We will identify P(x) as 
the fluid pressure. f2j is defined by 

= m A V in A , (18) 
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where m A is an auxiliary null vector such that m A n A = 
1. Using ([J|, for the black brane we find leading order 

fli = 2iTT v t . (19) 

We assume that the dynamics of the horizon geome- 
try perturbations are governed by the Einstein equations. 
The black brane is a solution to the Einstein equations 
with negative cosmological constant 

R AB + {d + l)9AB=Tf§ tt . (20) 

The Ricci tensor contraction R AB £ A e B can be expressed 
completely in terms of the horizon metric and extrinsic 
curvature [l2|. Consider first the component along n, 
which is the contraction with n A n B . This is the null fo- 
cusing equation, which after imposing the Einstein equa- 
tion takes the form 

-n A V A e+n(x)e-e 2 /d~a AB <J AB = T2 l B U n A n B . (21) 

There is no contribution from the cosmological constant 
term proportional to the metric due to |(2J). In the black 
brane case T™g U = 0. Plugging our previous results for 
the expansion, shear, and surface gravity into (|2ip . we 
find at leading order the incompressibility condition 

d iVi = . (22) 

Now consider the components of the Ricci tensor con- 
traction transverse to n. By definition, there is again no 
contribution from the cosmological constant term. From 
the remaining terms one finds 

C n Cii + d iK {x) - Dj4 + ^di6 = -n A e B T^ tt , (23) 
where 

C n a t = (d t + 9)VLi + r'D/.l, + VL 3 DiV j . (24) 

When T™g U — we find that the leading order terms 
are at order e 3 and give the NS equation (Q} without a 
force term and with a kinematic viscosity v = (47rT ) -1 . 
In cases where there is a non-zero matter stress tensor it 
will act as a forcing term in the NS equations, the force 
being 

fi = T£g*n B e£ . (25) 

For instance, adding a dilaton cj) to gravity, results in 

fi = V A 4>V B 4>n B e A /2 . (26) 

In the following we will consider the NS equations with- 
out a forcing term. From the NS equations one can derive 
the energy balance equation 

J \d t v 2 d d x = - J vd i v j d i v j d d x , (27) 

that relates the rate of change of the fluid energy to minus 
the energy dissipation per unit time due to fluid friction. 



To interpret this equation in terms of the horizon ge- 
ometry, consider the focusing equation (|21[) expanded to 
order e 4 with Tjftf * = 0, 

2ttT 6 = a AB a AB . (28) 

The expansion of the horizon is defined as the fractional 
rate of change in the cross-sectional area along the hori- 
zon generators 

6 = C n liiVh, (29) 

where C n is the Lie derivative and h is the determinant 
of htj . Using the incompressibility condition, (fl~3"|) for the 
shear, and ^ for the metric and then integrating l[28p 
over a horizon cross-section one has 

d t A = v(nT ) d J diVjdV d d x , (30) 

where A is the total horizon area. Imposing the energy 
balance law (|27[) . we find that 

d t (A/Aa) = - J d t v 2 /2d d x, (31) 

where Aq is the zeroth order area density (7rTo) rf - Thus, 
as the kinetic energy of the fluid on the boundary de- 
creases in time due to viscous dissipation, the horizon 
area grows. This is consistent with the classical area in- 
crease theorem of General Relativity [lj|. 

The derivation of the NS equations required knowledge 
of the horizon embedding and employed a local analy- 
sis near this horizon. Unlike [7], there was no need to 
know the asymptotic structure of the full bulk space- 
time. Therefore the results will apply to perturbations 
of a stationary non-singular null hypersurface, as long 
as there is a separation between the characteristic scale 
L of the perturbations and some intrinsic microscopic 
scale. The non-singularity requirement was used when 
contracting the Einstein equations in order to obtain the 
membrane equations. To see why a separation of scales 
is crucial, consider for example black holes in asymptot- 
ically flat spaces with horizons of a spherical topology 
(e.g. Schwarzschild). By dimensional analysis the corre- 
lation length of a fluid will scale as l c ~ Tq , where T 
is the Hawking temperature. In the asymptotically flat 
cases T^ 1 ~ ro, where Tq the horizon radius. Since the 
horizon is now also compact, the L can be no greater 
than ~ tq. Thus the dimensionless Knudsen number 
Kn = l c /L is of order unity, implying that the derivative 
expansion used above is not valid and that hydrodynam- 
ics is not the appropriate effective description. Moreover, 
in the general NS-like equation (f2"3")) found in the Mem- 
brane Paradigm approach, the term did does not vanish 
and leads to the a ssig nment of an unphysical negative 
"bulk viscosity" d, [If. 

To show that our analysis is applicable to cases other 
than the black branes in AdS, consider the example of 
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Rindler space associated with accelerated observers in 
(d + 2)-dimensional flat Minkowski space-time 

d 

ds 2 = -K 2 £ 2 dT 2 +df + J2 dxldx * > ( 32 ) 

i=l 

where k is a constant. To the uniformly accelerated ob- 
server with worldline £ = const., the surface £ = is a 
causal horizon, with intrinsic metric h%j — 5ij, that pre- 
vents him from an access to the entire space-time. The 
constant k can be identified with a temperature in the 
following way. Unruh [l|| showed that accelerated ob- 
servers feel the quantum vacuum to be a thermal state 
at temperature T — a/2n, where a the observer's proper 
acceleration. This Tolman local temperature can be ex- 
pressed as T — k/2itx, where x = V~9tt = k£ is the 
redshift factor. We define k — 2ttTq as the location inde- 
pendent temperature of the system. 

To perturb this horizon we allow, as in the black brane 
example, a slowly varying (L 3> Tq -1 ) fluid velocity n l = 
v % of ~ e. In EF coordinates the Rindler metric to O(e) 
is 

d 

ds 2 = -2irT Q rdt 2 + dtdr + V" dx i dx i 

i=l 

-Vidx*dr - AirT rdx l dt. (33) 

We assume the perturbations are governed by the vac- 
uum Einstein equations with zero cosmological constant, 
Rab — 0. Parameterizing k(x) as in (|17|) and ap- 
plying the membrane analysis to the Rindler horizon 
(r = 0) also shows that its dynamics is determined by 
the incompressible NS equations with kinematic viscos- 
ity (47rXb) _1 - In this case though the result cannot be 
understood as mirroring the hydrodynamics of a field the- 
ory fluid living on an asymptotic boundary of space-time. 
Here the fluid system could be interpreted as the vacuum 
thermal state (thermal gas of Rindler particles) perceived 
by accelerated observers, with the horizon acting as the 
holographic screen [17] . 

Consider now some of the implications of our results. 
First, exact solutions of the NS equations such as vortices 
and others [H.[l8| are mapped immediately into geometri- 
cal horizons. Second, finite time singularities, where fluid 
velocity gradients diverge, are mapped into the naked 
curvature singularities in the gravity description. Thus, 
the extremely important problem of singularities in the 
NS equations in three space dimensions appears linked 
to a cosmic censorship principle in gravity, a fact that 
may be of use in the future research of this long-standing 
problem. Third, in the case of turbulence, while an in- 
stantaneous picture of a flow is analytically unavailable, 
the spatially averaged correlators seem to have a well- 
defined analytic form. There is experimental and nu- 
merical evidence that in the "inertial" range of distance 



scales r much larger than characteristic length of viscous 
dissipation I but much smaller than flow excitation scale 
L, the flows exhibit a universal behavior, e.g. the space- 
averaged equal-time correlators of velocity differences in 
the inertial range are characterized by critical exponents. 
For instance, the longitudinal n-point functions scale as 

5 B (r) = <[(v(r)-v(0)).r/r] n )~^- . (34) 

The 1941 exact scaling result of Kolmogorov £3 = 1 
agrees well with the experimental data. In the geo- 
metrical picture, S n (r) correspond to the space-averaged 
equal-time correlators of differences of normals to the 
horizon. The fluid picture implies that in the evolution 
of the horizon surface there is a range where its area is 
preserved and it exhibits multi-fractal properties charac- 
terized by the exponents £„. 

Multi-fractality of the horizon surface may provide a 
dynamical basis for the so-called multifractal model of 
turbulence. This model, which appears to be the most 
successful phenomenological model of turbulence, ex- 
presses the anomalous exponents of turbulence in terms 
of the spectrum of fractal dimensions [3, [2l| . The model 
has been proposed originally on the basis of presumed 
singularities of the dynamical equations [l9| . However, 
it has been recognized that this assumption may not be 
correct @, [2(| ■ This leaves the open question of whether 
there is a singular structure underlying turbulence. A 
possible natural answer within the presented geometrical 
picture is that the horizon itself exhibits a multi-fractal 
structure in a wide range of scales. For instance, within 
the approximation corresponding to the Euler equation 
in hydrodynamics (infinite To), the dynamics may make 
the horizon surface indefinitely rugged, while preserving 
its area. This could correspond to the assumption of sin- 
gularities of the Euler equation. Thus, a central question 
is the horizon surface structure, which requires an under- 
standing of the general features of the surface fracturing. 
The geometrical picture is likely to be a useful tool also 
in the study of the relation between the velocity scaling 
in space and time, where in the latter one follows the La- 
grangian trajectory of a fluid particle in a turbulent flow 
and considers moments of velocity difference at time t 
1]. 
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